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Capillary instability and the resulting dynamics in an immiscible two-component Bose-Einstein 
condensate are investigated using the mean-field and Bogoliubov analyses. A long, cylindrical con- 
densate surrounded by the other component is dynamically unstable against breakup into droplets 
due to the interfacial tension arising from the quantum pressure and interactions. A heteronuclear 
system confined in a cigar-shaped trap is proposed for realizing this phenomenon experimentally. 
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I. INTRODUCTION 

A fluid cylinder will be unstable against breakup into 
droplets if its length exceeds its circumference. A fa- 
miliar example is a water jet projected into air, which 
breaks up into many water droplets. Savart [l| first in- 
vestigated this kind of instability and Plateau [2^ sub- 
sequently developed an experimental technique. Lord 
Rayleigh [3] used linear stability analysis to investigate 
this phenomenon. This instability is known as capillary 
instability \a\ (or sometimes as Plateau-Rayleigh insta- 
bility). The present paper proposes a method for observ- 
ing capillary instability in quantum degenerate gases. 

Classical fluid instabilities have recently been stud- 
ied theoretically in quantum fluids with renewed inter- 
est. Various fluid instabilities are predicted to be ob- 
servable in two-component Bose-Einstein condensates 
(BECs) including the Rayleigh-Taylor instability ^, |6|] , 
the Kelvin-Helmholtz instability [7, 8], the Richtmyer- 
Meshkov instability [9], the counter-superflow instabil- 
ity [lol, [Hi, and the Rosensweig instability [12]. Of 
these, the counter-superflow instability has recently been 
realized experimentally [l3|. The Benard-von Karman 
vortex street, which is predicted to occur in single- 
component Q and two-component BECs [l5|, is also 
caused by fluid instability. 

In the present paper, we investigate capillary insta- 
bility and the resulting dynamics in a two-component 
BEC. Capillary instability in a water jet originates from 
the surface tension of water that results from the attrac- 
tive interaction between water molecules. By contrast, 
in the present study, we demonstrate that capillary in- 
stability emerges in an atomic gas with repulsive inter- 
and intra-atomic interactions. We consider an immiscible 
two-component BEC in which a cylinder of component 
1 is surrounded by component 2. The interfacial tension 
between the two components [l6|, induces capillary in- 
stability and component 1 breaks into droplets that form 
a matter-wave soliton train. This system thus exhibits 
capillary instability and is novel in that the instability 
is only due to repulsive interactions. We also demon- 
strate that dark-bright soliton and skyrmion trains can 
be generated in this system. 

The present paper is organized as follows. Section Ull 



numerically demonstrates capillary instability and sub- 
sequent dynamics in an ideal system. Section IIIII ex- 
amines the stability of the system using Bogoliubov and 
linear stability analyses. Section |IVl proposes a realis- 
tic system for observing capillary instability in a trapped 
two-component BEC and numerically demonstrates vari- 
ous dynamics. Section |Vl presents the conclusions of this 
study. 



II. DYNAMICS IN AN IDEAL SYSTEM 



We consider a two-component BEC of atoms with mass 
rrij in an external potential Vj , where the subscript j = 1 , 
2 indicates component 1 or 2. In zero-temperature mean- 
field theory, the macroscopic wave functions ijji and ?/^2 
obey the Gross-Pitaevskii (GP) equation given by 
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where the interaction parameters are defined as 
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with ajji being the 5-wave scattering length between the 
atoms in components j and 2' {j^j' = l?^). We assume 
that the interaction parameters satisfy 



^11^22 < 912^ 



(3) 



for which the two components are immiscible [18|. 

In this section, we consider an ideal system with Vi = 
V2 = 0, mi = m2 = m, and gu = ^22 = 9- Normalizing 
the length and time in Eq. H]) by ^ = %/ (mgnY^'^ and 
^/'Ug, where n is the atomic density |'02p far from the in- 
terface and '^s = {gn/mY^'^ is the sound velocity, we find 
that the intracomponent interaction coefficients are nor- 
malized to unity and the relevant parameter is only gx^l g- 
The initial state is a stationary state ^ j (r) in which com- 
ponent 1 is a cylinder centered about the z axis and is 
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FIG. 1: (Color online) Capillary instability and resulting dy- 
namics of a two- component BEC for gi2/g =1.1. The three- 
dimensional objects depict the isodensity surface of compo- 
nent 1. The lower panels in (a)-(c) show the density distri- 
butions at the cross section of x = (the field of view is 
600^ X 100^). The number of component 1 atoms in an initial 
cylinder with length ^ is A^i^ = 400. 



surrounded by component 2. The state ^j(r) depends 
only on r = (x^ + ^^)^^^ and has cyhndrical symmetry 
and translation symmetry along z. Numerically, this ini- 
tial state is prepared by imaginary time propagation of 
the GP equation ([1]). The time evolution of the system 
is obtained by numerically solving the GP equation ^ 
by the pseudospectral method. 

Figure [1] shows the time evolution of the system for 
912/ g = l-l, where the number of component 1 atoms in 
an initial cylinder with length ^, 
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is 400. We add a small white noise to the initial state to 
break the exact numerical symmetry. At i = tv^/^ = 
3600, "varicose" [4] modulation with a wavelength 
99^ ^ 9.OR0 occurs on the cylindrical interface [Fig. [1] 
(b)], where Rq is the initial radius of the cylinder of com- 
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FIG. 2: (Color online) Time evolution of the system with a 
periodic initial perturbation given by Eq. ([5]). The parameters 
are the same as those in Fig.[Tl The three-dimensional objects 
indicate the isodensity surfaces of component 1 and the lower 
panels show the density distributions at the cross section 
of X = (the field of view is 600^ x 1000- 



ponent 1. The varicosities then split into droplets [Fig.[T] 
(c)], which undergo quadrupole oscillations [Figs. [T] (d) 
and [H (e)] due to the kinetic energy converted from the 
interfacial tension energy. The droplets subsequently ag- 
gregate to minimize the interfacial tension energy. 

Figure [2] shows the time evolution of the system, where 
a periodic initial perturbation is added to component 1 
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instead of the white noise. This wavelength 200^ is about 
twice the most unstable wavelength in the modulation in 
Fig. [1] (b). The sinusoidal seed develops into a highly 
nonlinear pattern [Fig. [2] (a)] and the thin ligaments sep- 
arate from the main droplets [Fig. [2] (b)], which deform 
into small drops [Fig. [2] (c)]. Such nonlinear behavior is 
very similar to that in classical liquid jets [19j, and the 
small drops are termed "satellite drops" . 



III. STABILITY ANALYSIS 

A. Bogoliubov analysis 

We study the stability of an ideal system such as that 
depicted in Fig. [1] (a) using Bogoliubov theory. We sep- 
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FIG. 3: (Color online) Radial density profiles of the sta- 
tionary states for (a) A^i^ = 400 and (b) iVi^ = 7854 with 
912/9 = 1.1. (c) Imaginary part of the excitation frequencies. 
The circles and triangles are obtained by numerical diago- 
nalization of Eq. (|9)) with the parameters used in (a) and 
(b) , respectively. The plotted quantities are normalized using 
Eq. ([To)) with i?o = for the circles and i?o = 50^ for the 
triangles. The solid curve indicates Eq. (|30p . 



Vi(r,^,t) = [vl>,-(r) + 5i^^{r,z,t)] e-'^^■*/^ (6) 



where fij is the chemical potential and the system is as- 
sumed to have cylindrical symmetry. Figures [3] (a) and [3] 
(b) show examples of the stationary state ^j(r). Substi- 
tuting Eq. (|6j) into the GP equation ([1]) and taking the 
first order of the small deviation we find (j ^ j') 
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where and ^2 are assumed to be real without loss of 
generality. Using the excitation mode of the form 



in Eq. ([7]), we obtain the Bogoliubov-de Gennes equation 
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Figure [3] (c) shows the imaginary part of the Bogoli- 
ubov excitation frequency obtained by numerical diago- 
nalization of Eq. (|9]). For comparison with an analytic 
result given in the next subsection, we plot the normal- 
ized quantity as 



Imcj = \l ^Imcj, 



(10) 



where a is the interfacial tension coefficient originating 
from the excess kinetic and interaction energies at the 
interface. For 912/9 = 1.1, we can use the expression for 
the interfacial tension coefficient obtained in Refs. [TgI, 
fl7l |. namely, 
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Figure [3] (c) reveals that the most unstable wavelength is 
:^ 9.3i?05 which agrees with the modulation in Fig.[T](b). 
In Fig. [3](c), we plot two cases for the parameters corre- 
sponding to Fig. [3] (a) (circles) and Fig.[3](b) (triangles). 
We note that these plots fall on a universal curve using 
the scaling in Eq. (p^Oj) . 



B. Linear stability analysis 

To understand the Bogoliubov spectrum obtained in 
Fig. O we analyze the stability of the system using a 
similar method to Rayleigh's linear stability analysis [sl. 
We start from the mean- field Lagrangian given by 
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where 



Pj = mr-^ - ^iWii' - vm' - ^i^.f • (13) 

The functional derivative of the action J Ldt with respect 
to gives the GP equation ([1]). In this subsection, we 
assume that the two components are strongly phase sepa- 
rated and that the interface thickness is negligibly small. 
We also assume that the system has cylindrical symme- 
try and that the interface is located at r = R{z^t). We 
then approximate Eq. ([12]) as 



L = 27T 



j dz rdrPi + j rdrPs^ - aS, (14) 



where S is the area of the interface. Differentiating 
Eq ([Hj) with respect to we obtain [20] 

Pi{R,z,t) - P2{R^z,t) = a (^i- + i-^ , (15) 

where Ri and R2 are the principal radii of the interface 
curvature. Equation ([15]) corresponds to Laplace's for- 
mula in fluid mechanics. 

We write the wave function as 



^Pj{r,z,t) = ^nj{r,z,t)e"f^'-^'^''\ (16) 

where rij and are real functions, ni = for r > R, and 
722 = for r < R' We separate the density and phase as 



We assume the sinusoidal forms of the small excitation 

as 

R{z,t) = Ro^esm{kz-ujt), (23a) 
Vjr{r^z^t) = Vjr{r) cos{kz — ut), (23b) 
Vjz{r^z^t) = Vjz{r) sm{kz — ujt)^ (23c) 
Suj = Shj{r) sm{kz — ujt). (23d) 

Using Eqs. ([22]) and ([23]), Eq. ([IS]) becomes 

giinioShi{r) - g22n2oSh2{r) = ^{k'^Rl - 1), (24) 

where we have used the equilibrium condition (^nn^o — 
^22n^o)/2 = a/Ro. Substituting Eqs. (|23bl) - (l2Ml) into 
Eqs. ([IH]) and ([21]) yields 

r^vjrir) + rv'j^{r) - (1 + k\^)vjr{r) = 0. (25) 

The solutions of this differential equation are 

vir{r) = ci/i(/cr), (26a) 
V2r{r) = C2Ki{kr), (26b) 

where ci and C2 are constants and In and are modified 
Bessel functions of the first and second kinds. These 
functions are chosen such that components 1 and 2 do 
not diverge at the z axis and infinity, respectively. It 
follows from Eqs. ([21]) and ([26]) that 
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nj(r,z,t) = Hjo + Srijir, z,t), 
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4>j{r,z,t) = -i21^t + 5Mr,z,t), (17b) 



and substitute them into Eq. ([!]). Taking the first order 
of S<pj and gjjSrij^ we have 



V-Vj = 0, 

, dS(l)j 

n-g^^gjjSrij = 0, 
where the velocity is defined as 



Vj — VOOj. 



The gradient of Eq. ([T9|) gives 
Ov ■ 



(18) 
(19) 

(20) 
(21) 



Substituting Eq. ^ into Eq. and using Eq. p^ . 
we obtain the pressure, 



(22) 



The kinematic boundary condition at the interface is 
dR/dt = vir{R) = V2r{R)^ which gives 



-cih{Ro) = -C2Ki{Ro). 



(28) 



Using Eqs. ([24]), ([27]), and ([28]), we obtain the disper- 
sion relation of the excitation as 



kRoiPRl - 1) 



(29) 



This dispersion relation has the same form as that of 
classical inviscid incompressible fluids [2l| . The quantum 
mechanical pressure is included in the interfacial tension 
coefficient a and Eq. (|29|) does not contain any explicit 
quantum correction term. For < kRo < 1, the right- 
hand side of Eq. ([29]) is negative and the frequency u 
is pure imaginary. The mode with a wavelength larger 
than 27tRo is therefore dynamically unstable. The most 
unstable wave number is given by kRo ^0.68 and the 
most unstable wavelength is ~ 9.2i?o, which is in good 
agreement with the modulation wavelength in Fig.[T](b). 
The most unstable wave number depends only on Rq and 
is independent of the interaction parameters. The inter- 
action, which is included in a, only affects the growth 
rate of the unstable modes. 
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for < kRo < 1, which is in excellent agreement with 
the numerically obtained Bogoliubov spectrum normal- 
ized by Eq. (p!Q|) . This confirms that the modulation in- 
stability shown in Fig. [T] (b) is the capillary instability. 
Figure [3] shows that the dispersion relation ([29|) is accu- 
rate even when the radius of the cylinder of component 
1 is of the same order as the interface thickness. 



IV. DYNAMICS IN A TRAPPED SYSTEM 




We propose a realistic experimental situation to ob- 
serve the capillary instability in a trapped BEC. We con- 
sider a heteronuclear condensate consisting of ^^K (com- 
ponent 1) and ^'^Rb (component 2) atoms in the stretched 
\F = 2, mi? = 2) states. [22]. The scattering lengths are 
ail = 65aB 0, (222 = 99aB |2J], and ai2 = 163aB pBj . 
which satisfy the phase separation condition in Eq. (jSj. 
The two components are confined in axisymmetric har- 
monic potentials, 

V, = ^{o.lr'+oj%z% (31) 

produced by laser beams. Since the electronic excitation 
frequencies of the two atomic species are separated from 
each other, Vi and V2 can be chosen independently by 
using laser beams with different frequencies. 

The initial state is prepared for trap frequencies sat- 
isfying uji± ^ UJ2± and ujj± ^ ujjz. Since the radial 
confinement of component 1 is much tighter than that 
of component 2, the ground state of component 1 is ex- 
tremely narrow and it is surrounded by component 2 in 
the radial direction. The tight radial confinement pre- 
vents component 1 from breaking up into droplets for the 
same reason as why a classical liquid in a pipe never ex- 
hibits the capillary instability. We then add small white 
noise to the initial state and reduce uji± suddenly at 
t = 0. This produces an unstable initial state similar 
to the state in Fig. [T] (a), and capillary instability is ob- 
tained. 

Figure [4] demonstrates the capillary instability and re- 
sulting dynamics of a trapped two-component BEC ob- 
tained by numerically solving the 3D CP equation. The 
ratio of the number of atoms in component 1 to the total 
number of atoms is Ni/{Ni + N2) = 0.05. At t = 0, the 
radial confinement of component 1 is relaxed and mod- 
ulation grows due to the capillary instability. Varicose 
modulation occurs in component 1 [Fig. |4] (b)], which 
breaks up into droplets [Fig. |4] (c)]. The axisymmetry 
of the system is preserved within this time scale. We 
have thus shown that the capillary instability can be ob- 
served in a trapped two-component BEC. The structure 



FIG. 4: (Color online) Time evolution of a ^^K-^'^Rb conden- 
sate confined in a cigar-shaped trap. The initial state is the 
ground state, in which Ni = 5 x 10^ ^^K atoms (component 
1) are confined in Vi with {(jJi±,(jJiz)/27v = (582, 11.6) Hz and 
N2 = 9.5 X 10^ ^^Rb atoms (component 2) are confined in V2 
with {(j02±,oj2z)/27v = (89.4,7.2) Hz. At t = 0, cji±/(27r) is 
changed from 582 Hz to 145 Hz. The panels in (a)-(c) show 
the density distributions iV^jl^ at the cross section of x = 0, 
where the unit of the density is 10^° m~^ and the field of 
view is 220 x 14 /xm. The bottom figure shows the isodensity 
surface of component 1 at t = 21 ms. 



in Fig. m (c) is similar to a matter- wave soliton train. A 
matter-wave soliton and soliton train have so far been 
produced by a single-component BEC with attractive in- 
teractions [26, 27]. By contrast, the train of matter- wave 
droplets in Fig. S] (c) is produced only by repulsive inter- 
actions. 

Figure [5] shows the dynamics for Ni/{Ni + A^2) = 0.2. 
Since the radius of the initial cylinder of component 1 is 
larger than that in Fig. SI the modulation wavelength for 
the capillary instability is larger. The size of each droplet 
produced by the capillary instability is thus larger than 
that in Fig. |4] and the droplets cut component 2 along 
the z axis [Fig. [5] (c)]. The two-component structure 
in Fig. [5] (c) is similar to dark-bright solitons ^28|-(30|, 
where the density dips in component 2 and the droplets 
of component 1 correspond to dark and bright solitons, 
respectively. In fact, the phase 02 of component 2 jumps 
at the density dips [dashed square in Fig.[5](c)], which is 
required to stabilize dark-bright solitons. 

Figure [6] demonstrates generation of a skyrmion train 
using the capillary instability. The initial state is a sta- 
tionary state in which component 1 has a singly quan- 
tized vortex line along the z axis [Fig. [6] (a)]. Such a 
state is generated by phase imprinting using, for example, 
a stimulated Raman transition with Laguerre-Gaussian 
beams At t = 0, the radial confinement of compo- 
nent 1 is reduced and a force in the z direction is exerted 
on component 1; this can be achieved, for example, by 
shifting the trap center in the z direction. Figures [6] (b) 
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FIG. 5: (Color online) Time evolution for Ni 
and A^2 = 8 X 10^. The trap frequencies are (cji^, 
(436,35) and {uj2±,uj2z)/27t = (89.4,7.2). At t = 0, uji±/27v 
is changed to 145 Hz. The panels show the column density 
dj — f \ipj\'^dx and the phase 0^ = arg'0j(x = 0). The unit 
of the column density is 10^^ m~^ and the field of view is 
170 X 13.6 /im. The dashed square in (c) shows an example 
of a position at which the phase jumps. 



and[6](c) show the state at t = 9 ms. The isodensity sur- 
face in Fig.[6](c) indicates that the initial hollow cylinder 
of component 1 breaks up into toruses due to the capil- 
lary instability. The initial vortex remains in component 
1, as shown in Fig. [6] (d). This breakup is not due to 
the Kelvin-Helmholtz instability since it occurs in a very 
similar manner when there is no force in the z direction. 
The vortex rings are generated in component 2 [Fig. [6] 
(b)] at locations where there are toruses of component. 
This vortex ring generation is due to the acceleration of 
the toruses of component 1 in the z direction [15]. The 
structures produced in Fig. [6] are thus the skyrmions in a 
two-component BEC [32]. However, the skyrmion train 
in Fig. [6] is unstable against axisymmetry breaking. Be- 
cause of the vortex quantization, the behavior in Fig. [6] 
is quite different from that in an axially rotating jet of a 
classical fluid 1331. 



V. CONCLUSIONS 

We have investigated the capillary instability and re- 
sulting dynamics in a two-component BEC. We first 
demonstrated the dynamics in an ideal system in Sec. |TT1 
We have shown that modulation occurs on a cylinder 
of component 1 due to the capillary instability and the 
cylinder breaks up into droplets (Fig. [T]), as in a classi- 
cal fluid jet. Formation of satellite drops was also ob- 











z 

(b) t = 9 ms 


'^2' 1 






















ft § 1 1 i 











(c) 



ooooo 



A, 




FIG. 6: (Color online) Creation of a skyrmion train, (a) 
Density distribution of the initial state at the cross section 
of X = 0. Component 1 has a singly quantized vortex along 
the z axis. The initial trap frequencies for component 1 are 
(cJi±,cJiz)/27r = (796,14.5) Hz. At t = 0, uJi± is changed 
to 145 Hz and the trap center for component 1 is shifted to 
z = 65 /im. The other parameters are the same as those in 
Fig. m The field of view is 150 x 14 /xm. (b) Density profile 
of component 2 at the cross section of x = at t = 9 ms. 
The density and phase in the dashed square are magnified 
in the lower panels. The arrows indicate the directions of 
circulations of the quantized vortices, (c) Isodensity surface 
of component 1 at t = 9 ms. The lower panels show the 
density and phase profiles of component 1 at the cross section 
of z indicated by the arrow. The unit of the density in (a)-(c) 
is 10^° m"^ 



served (Fig. [2]). In Sec. ITTH we performed Bogoliubov 
analysis and numerically obtained a dynamically unsta- 
ble spectrum, which is in good agreement with that ob- 
tained by Rayleigh's linear stability analysis (Fig. [3j). In 
Sec. HVl we proposed realistic trapped systems and nu- 
merically demonstrated the dynamics caused by the cap- 
illary instability. We have shown that the capillary insta- 
bility can be observed in a heteronuclear two-component 
BEC confined in a cigar-shaped harmonic trap, where 
the radial confinement of inner component is controlled 
(Fig. 3]). We have also shown that dark-bright soliton 
and skyrmion trains can be generated in this system, 
which are quantum mechanical objects and have no coun- 
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terparts in capillary instability in classical fluids. 
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